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We investigate a chain of capacitively coupled Josephson 
junctions in the regime where the charging energy dominates 
over the Josephson coupling, exploiting the analogy between 
this system and a multi-dimensional crystal. We find that 
the current-voltage characteristic of the current-driven chain 
has a staircase shape, beginning with an (insulating) non- 
zero voltage plateau at small currents. This behavior differs 
qualitatively from that of a single junction, which should show 
Bloch oscillations with vanishing dc voltage. The simplest 
system where this effect can be observed consists of three 
grains connected by two junctions. The theory explains the 
results of recent experiments on Josephson junction arrays. 

PACS numbers: 73.23.-b,74.50.+r 



Recent experiments carried out on two-dimensional 
Josephson junction arrays in the quantum regime - i.e. 
where the charging energy exceeding the Josephson cou- 
pling introduces quantum dynamics - showed remark- 
able steps of the current-voltage characteristics. Two 
voltage steps could be clearly seen, with voltage value 
given by the superconducting gap 2A/e. Further struc- 
ture at higher voltages is washed out. The purpose of 
the present article is to show that these steps can be 
naturally explained in the framework of the Bloch os- 
cillation description of quantum mechanical Josephson 
junctions. Moreover, we argue that the physics of a 
chain differs qualitatively from that known from single 
Josephson junctions ^ n contrast to a current-biased 
single-junction, a chain of ID Josephson junctions has no 
zero- voltage state, i.e. the average voltage is greater than 
2A/e for any current. The only stable state of the system 
has a voltage equal to 27V A/e where N is the number of 
junctions. However, for low currents this stable state is 
reached only on astronomical time scales, 

t z = 2er 1 exp(Iz/I), 

while the time required for the voltage 2A/e to establish 
is 



2eA/IE c <.t z . 



Here Ec is the charging energy, and Iz is defined be- 
low. Thus, in realistic experiments metastable I—V char- 
acteristic are observed, with voltage steps starting from 
the value 2A/e. The time required to achieve the stable 
voltage depends strongly on the external current. These 



results are a consequence of a lowered symmetry of the 
system, and the effect shows up already in the system 
shown on Fig. 1. We conjecture that the physics of 2D 
Josephson junction arrays is similar, which provides an 
explanation of the experiments 0. 

First we recollect the results of a single Josephson junc- 
tion with a capacitance C and Josephson coupling energy 
Ej. In the absence of dissipation and for energies below 
the superconducting energy gap A it can be described by 
the Hamiltonian 0,01 



H = — -Ejcost 



= -2ie- 



d 



(1) 



where Q and <p are the charge and the superconducting 
phase difference operators, respectively. The Hamilto- 
nian (|l) is equivalent to that of a quantum particle with 
a coordinate ^ in a periodic potential. If the junction is 
driven by an external current / (which is an analogue of 
an external force for a quantum particle), the junction 
charge increases with time, Q = It, until it reaches the 
value Q — e. Then the system has two options: either a 
Cooper pair tunnels across the junction and the junction 
charge is reset from Q — e to Q — —e, or if this does 
not happen the charge increases further with time. The 
former option is analogue to a Bragg reflection when a 
quantum particle reaches the edge of the Brillouin zone 
and is reflected, thus staying in the lowest energy band. 
The latter option is the equivalent of Zcncr tunneling 
across the band gap to a higher band Q. In the limit 
Ej -C Ec — e 2 /2C (which will be discussed in this pa- 
per) the probability of such Zener tunneling event is H 



A = exp(-J z /J), I z = neEj/8E c 



(2) 



The question is which of the two options are realized and 
how they manifest themselves in an experiment. The 
answer depends on the current: 

Low currents: In this regime the average time tz re- 
quired for a Zener tunneling event to occur exceeds the 
time of the experiment, tz = 2e(/A)~ 1 ^> t e . This means 
that during the whole experiment the system stays in the 
lowest energy band, the charge (and the voltage) across 
the junction oscillates with the period 2e/I around zero. 
This is the metastable regime of Bloch oscillations . 

High currents: In this regime, t e 3> tz, Zener tunnel- 
ing event || takes place during the experiment. In the 
absence of dissipation the system - after it jumps to the 
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second energy band - will continue to rise in energy |7]] 
until its energy reaches the value 2A + Ec (the average 
voltage across the junction at this point is Vq f» 2A/e 
A further increase is not possible because of dis- 
sipation due to single electron tunneling that sets in as 
soon as the voltage exceeds the value 2A/e. Then the 
voltage will oscillate around the stable value 2A/e ||. 

Thus the difference between the two above regimes 
is quite clear: the first one corresponds to metastable 
zero-voltage state of the junction whereas in the sec- 
ond regime a finite stable voltage V « 2A/e is mea- 
sured. The crossover between these two regimes occurs 
at / « I z /\n(It e /2e). 

Experimentally it is more convenient to study charging 
effects in systems containing several (or many) Josephson 
junctions. Hence one can ask whether the simple physical 
picture discussed above remains valid also in the case 
of two or more coupled Josephson junctions. Intuitively 
one could expect that for sufficiently low currents flowing 
through the chain of Josephson junctions the system will 
stay in the zero voltage state V ~ similarly to the case 
of a single junction. 

In this paper we argue that the behavior of already two 
capacitively coupled Josephson junctions is qualitatively 
different: a finite voltage drop V ^ Vo = 2A/e across the 
system will be measured no matter how low the external 
current is. This conclusion remains valid for an arbitrary 
number of junctions in the chain N > 1. To understand 
this result we use the analogy between our system and 
an electron in a multidimensional crystal. If the self- 
capacitance Co of the superconducting islands is non- 
zero the symmetry of the corresponding Brillouin zone is 
lowered in a way which prevents the system from staying 
in the lowest energy band for any non-zero current I. 
Zener tunneling in at least one of the junctions occurs 
immediately after the current is applied and the system 
switches to a finite voltage state. The resulting current- 
voltage characteristics of the Josephson junction chain 
has a form of a staircase. Our results can explain the low 
voltage behavior of two-dimensional Josephson junction 
arrays recently reported in Ref. Q . 
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(b) 

FIG. 1. (a) Two capacitively coupled Josephson junctions; 
(b) An experimental realization of the same system. 

For the sake of simplicity we first consider the system of 
three superconducting grains with phases 0j, i = 1,2,3, 
connected by two Josephson junctions (Fig. la). The 



relevant variables for description of this system are linear 
combinations of the phases (cf. e.g. ||) 

01 = 01 ~ 025 02 = 02 - 03; #3 = (01 + 02 + 03)/3. 

We can also introduce the junction charge operators 
Qi = —2ied/d8i, i — 1,2, and the "center-of-mass 
charge" Q 3 = —2ied/d9 3 . Then the Hamiltonian reads 
as 

6 = W + \ E Qiic-^aQi - Ej ]T cos ft, (3) 

ij i 

and the effective 2x2 capacitance matrix C has elements 
Cu = C 22 = C + 2C /3; C 12 = C 21 = C /3. The 
center-of-mass dynamics is independent of that of both 
junctions. The latter is essentially the motion of an elec- 
tron in a 2D periodic potential. In the limiting case 
Ej <C e 2 /max{Co,C} (nearly-free-electron model) the 
Josephson term in the Hamiltonian (||) can be consid- 
ered as a perturbation. We assume furthermore that 
2A > e 2 / max{Co, C}. The eigenfunctions of the un- 
perturbed Hamiltonian are "plain waves" , 

^o(0 ) = exp(iQ0 /2e + iQ 3 6 3 /2e). 

Here we have introduced two-dimensional vectors 9 = 
(0i,0a) andQ = (Qi,Q 2 ). 

Now we can clarify the meaning of the charges Qi 
and Q 2 . Introducing the grain charge operators qt = 
—2ied I d<pi, i = 1,2, 3, with the eigenvalues qi, we obtain 

qi = Ql+Q3/3; q 2 = -Q1+Q2+Q3/3; q 3 = -Q2+Q3A 

The charge Q 3 is just a total charge of the chain. It is 
conserved, and from now on we put it equal to zero jjj. In 
the current-bias regime the charge q\ grows linearly with 
time, except for the jumps by 2e due to Cooper pair 
tunneling (CPT) through junction 1 (see below). The 
charge q 3 decreases linearly with time, except for CPT 
through junction 2, while the charge q 2 keeps constant 
except for CPT through the either junction. In other 
words, the system can move only along lines Q\ — Q 2 = 
2en, with n being an integer, and CPT processes are 
responsible for the jumps of the system between these 
lines. The time evolution of the system starts from the 
line Qi — Q 2 — It. 

The eigenvalues of the unperturbed Hamiltonian are 

So(Q) = ^Q i (C- 1 ) ij Q i . (4) 

ij 

The Josephson coupling plays a role only in the vicin- 
ity of "critical surfaces", -E'o(Q) = Eo(Q — K), with 
K = (2pe, 2eq) (p and q are integer numbers) being the 
reciprocal lattice vectors. In lowest (second) order in Ej 
only vectors K with p = ±1, q = and p — 0, q — ±1 are 
essential, and the Josephson coupling creates gaps equal 
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to Ej on the critical surfaces. For the spectrum (|4|) these 
are given by 

1) Qi - aQ a = ±e; 2) Q 2 - aQ x = ±e, 

with a = (C Q /3)(C + 2C0/3)- 1 . Note that < a < 1/2; 
the limits a — and a = 1/2 correspond to the cases 
Co = and C = 0, respectively. Thus, the critical sur- 
faces are just two pairs of parallel straight lines. The 
romb formed by these lines is the first Brillouin zone. 
The Bragg reflection processes are possible on the crit- 
ical surfaces, i.e. the system can jump by a vector K. 
The energy is periodic in Q-plane with two orthogonal 
periods equal to 2e. It will be convenient for us to use the 
extended band picture. The voltage across the junction 
i is given by Vi — dE/dQi. 
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FIG. 2. The motion of the current-biased system in the 
plane (Qi,Q2); Co = 0. The critical surfaces are shown by 
boldface lines; the first Brillouin zone is a square ABCD. In 
the weak current regime the system moves along the line AC. 
Possible trajectories for the cases of strong currents, A ~ 1, 
are shown by dotted lines. 

For Co = the Hamiltonian ([|) is just a sum of two 
Hamiltonians ([!]) of isolated junctions, and the behavior 
of the chain is trivial. The critical surfaces are Qi = ±e, 
i = 1,2 (Fig. 2). For low currents a possible time evo- 
lution (with large probability (1 — A) 2 ) is simultaneous 
Bloch oscillations in two junctions ]lC|| , i.e. Bragg reflec- 
tion from point A= (e, e) to point C= (— e, — e). The 
mean voltage is zero for this process. Another possibil- 
ity (with probability A(l — A)) is the reflection to the 
point B= (— e, e), with subsequent motion along the line 



(Q, Q+2e), Q > — e. This means Zener tunneling injunc- 
tion 2, and Bloch oscillations in junction 1. Eventually 
the system reaches the critical surface in the point (e, 3e), 
where it is reflected to the point (— e, 3e) and so on. The 
voltage across junction 1 V\ = Qi/C oscillates with fre- 
quency I /2e, while the voltage across junction 2 grows 
linearly in time. The third possibility (with equal proba- 
bility A(l — A)) is an equivalent process, where the system 
is initially reflected to the point (e, — e). Finally, the last 
possibility (with probability A 2 ) is to proceed along the 
line (Q, Q) without any reflection. This corresponds to 
Zcncr tunneling in two junctions simultaneously, i.e. the 
voltage grows across both junctions. 
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FIG. 3. The same as Fig. 2, Co / 0. The first Brillouin 
zone is a romb ABCD. In the weak current regime the system 
moves along the dashed line BF. 

The behavior is quite different in the presence of a fi- 
nite self-capacitance Co- The system again starts from 
the point (0, 0) and reaches the critical surface in the 
point A= (e&eO, £ = (1 - 1 < € < 2 (Fig. 3), 

beyond the point (e, e). Consequently, it has only three 
possibilities of further motion. The first one (with proba- 
bility (1 — A 2 )/2, where A is given by the same expression 
(|) with E c = e 2 /2(C + C )) is to be reflected to the 
point G= (e(£ — 2),e£). It corresponds to Zener tun- 
neling in junction 2. The subsequent evolution of the 
system is voltage growth across junction 2 and Bloch os- 
cillations in junction 1. The equivalent possibility is to 
be reflected to the point H(= e£,e(£ — 2)). Finally, the 
last allowed process (with probability A 2 ) is to continue 
the motion along the line (It, It) without any reflection, 
i.e. Zener tunneling takes place in both junctions. Note 
that the simultaneous occurrence of Bloch oscillations in 
both junctions is impossible. This is a consequence of a 
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lowered symmetry: while the case Co = has the same 
symmetry group as a square, the introduction of Co low- 
ers the latter to that of a romb. 

Thus, the situation differs qualitatively from that of a 
single junction. Still two scenarios exist, low and high 
currents. The second one yields the mean voltage 2Vb, 
which one could expect. However, the low current sce- 
nario leads now to a finite mean voltage Vq , in contrast to 
the single-junction case. The time t s ~ 2eA/IEc <§C tz, 
required for this metastable value to establish, is very 
short. Hence at the beginning of experiment the mean 
value Vq is observed. The time of order tz, which can 
be extremely long, is required to obtain the stable value 
2V . 

if 



1 2 3 V/V 

FIG. 4. A sketch of voltage-current characteristic of the 
system of coupled junctions. All fine structure of oscillations 
is eliminated, and only the mean value V in units of Vq is 
shown. 

The above approach can be generalized to the case of 
the Josephson chain consisting of N + 1 grains and N 
junctions. The relevant variables are phases of junctions 
6i = 4>i — 4>i+i, 1 < i < N, and the center-of-mass co- 
ordinate 9 N+ i = (N + iy 1 ((f>i + <p 2 + ■ ■ ■ + 4>n+i), with 
4>i being the phases of the grains. The problem is equiv- 
alent to the motion of an electron in an iV-dimensional 
periodic potential; the chain is represented by a point 
Q in the space of junction charges Qi, i — 1, . . . , N. A 
charge-neutral system moves along the lines 

(Qi + 2eNi, Q 2 + 2en 2 , . . . , Qn + 2en N ), n k E Z. 

It starts from the coordinate origin and can jump be- 
tween the lines only as a result of "Bragg reflection" on 
the critical surfaces, given by (I = 1, . . . , N) 

^Qi^)* = ^(^U' (5) 

i 

where the N x N capacitance matrix is given by 



Cm = C lk = C8 kl + C v - - - ' ; k < I. (6) 

An analysis for arbitrary relation between capacitances 
C and Co would require the inversion of the capacitance 
matrix C. However, both cases C = and Cq <C C 
allow an analytical solution and show the same symmetry 
lowering causing a finite (metastable) voltage Vb across 
the chain. 

Thus, the current-voltage characteristics can be de- 
scribed as follows. The stable value of the mean volt- 
age for any value of current is NVq. However, the time 
required to achieve this stable voltage is astronomical 
(exponentially long) for low currents, and may exceed 
the duration of the experiment t e . Hence one observes 
the metastable value Vq. The main effect of an increasing 
current is a decrease of the time tz- For high currents one 
has tz -C t e and observes the stable voltage. There are 
some intermediate regimes between, and metastable volt- 
ages kVo, k = 2, 3, . . . , N — 1 can be measured (Fig. 4). 
This dependence is essentially a low-voltage that ob- 
served in the experiments . 

In conclusion, we predict that a current-biased chain 
of N capacitively coupled Josephson junctions always is 
in a finite voltage V > 2A/e state no matter how low the 
external current is [0. This conclusion differs qualita- 
tively from that reached for a single Josephson junction. 
It is based on simple symmetry arguments. While for 
vanishing on-site capacitance Co = the structure of the 
Brillouin zone for a Josephson "particle" corresponds to 
that for an N-dimcnsional cubic crystal, the crystal sym- 
metry is lowered for any non-zero Co- As a result Bloch 
oscillations in the lowest Brillouin zone cease to occur 
and the system "jumps" to higher bands due to Zener 
tunneling. This value of the voltage is metastable, but 
can be measured experimentally since the time scale re- 
quired to observe the stable value NVq is very long fol- 
low currents. We found a staircase- like shape of the I-V 
curve of the system and argue that our predictions can 
explain the observed low-voltage behavior of Josephson 
junction arrays pi. 
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